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Some past and ongoing explorations of the spectrum of QCD using Monte Carlo simu- 
lations on a space-time lattice are described. Glueball masses in the pure-gauge theory 
are reviewed, and the energies of gluonic excitations in the presence of a static quark- 
antiquark pair are discussed. Current efforts to compute the baryon spectrum using 
extended three-quark operators are also presented, emphasizing the need to use irre- 
ducible representations of the cubic point group to identify spin quantum numbers in 
the continuum limit. 

Keywords: Lattice QCD calculations; Lattice gauge theory; Glueballs; Baryons. 
PACS numbers: 12.38.Gc, ll.15.Ha, 12.39.Mk 

1. Introduction 

Spectroscopy is a powerful tool for distilling the key degrees of freedom in a given 
system. At the present time, the best way of extracting the spectrum of states 
from the QCD Lagrangian is Monte Carlo computer calculations using a space-time 
lattice. A spectrum determination requires the extraction of many excited-state 
energies, so a brief discussion on how excited-state energies can be determined 
from Monte Carlo estimates of correlation functions in Euclidean field theory is 
warranted. In this talk, several past explorations of the QCD spectrum are outlined, 
in particular, glueball masses in the pure-gauge theory and the energies of gluonic 
excitations of the so-called static quark-antiquark potential. Also, ongoing efforts by 
the Lattice Hadron Physis Collaboration (LHPC) to determine the baryon spectrum 
using extended three-quark operators are described. 

2. Extracting Excited-State Energies and Resonances 

The Monte Carlo method can be applied to obtain estimates of the path integrals 
which yield a Hermitian matrix of correlation functions Cy ;(f) = (0| Oi(t)Oj(0) |0), 
where Oj(0) creates the states of interest at time t = and Oi(t) annihilates 
such states at a later time t. The procedur e for extracting the lowest energies 
Eq, E\, E2, ■ . ■ from this matrix is well knowrP^. Let A n (t,io) denote the eigen- 
values of the hermitian matrix C(^o) -1 ^ 2 C(t) C(to) -1 ^ 2 , where to is some fixed 
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reference time (typically small) and the eigenvalues, also known as the principal 
correlation functions, are ordered such that Ao > Ai > • • • as t becomes large. Then 
one can show that 

Km A„(Mo) = e- £ " (t -' o) (l + 0(e- A " (t - to) )), (1) 
A n = min|.E fc -.E n |. (2) 

Determinations of the principal correlators X n (t,t ) for large temporal separations 
t yield the desired energies E n . Since statistical fluctuations grow with increasing 
t, it is crucial that contributions from higher-lying states be diminished so that 
the desired lowest-lying energies dominate the principal correlators well before the 
signal-to-noise ratio falls. Judiciously chosen quark-field and gluon-field smearings 
is one important ingredient for reducing couplings to the short wavelength modes 
of the theory. The use of large sets of extended operators is another key ingredient. 

Our Monte Carlo calculations are carried out in a finite-sized box with periodic 
boundary conditions. Given the finite volume and the discrete nature of the allowed 
momenta in the box, the masses and widths of resonances (unstable hadrons) cannot 
be calculated directly, but must be deduced fro m the d iscrete spectrum of finite- 
volume stationary states for a range of box Such applications require 
prohibitive computational resources, but our goal in the baryon project is to ob- 
tain a first exploratory scan of the spectrum, so simply obtaining the finite- volume 
spectrum for a few judiciously-chosen volumes should suffice for ferreting out the 
hadron resonances from the less interesting scattering states. 

3. Glueballs and Gluonic Excitations in Presence of Static QQ Pair 

Glueball masses in the pure-gauge theory, shown in Fig. ^ were the first exploration 
of the QCD spectrum^ in which I was involved. The spectrum can be fairly well 
described by a bag model description of the gluons, whereas string models seem to 
fare less well. Inclusion of light-quark effects is an ongoing challenge. 

The issue of string formation of the gauge-field in the presence of a static quark- 
antiquark pair was my next exploration . For quark-antiquark separations R greater 
than 2 fm, the spectrum of gluonic excitations agreed without exception with that 
expected from an effective string theory description of the gauge field (see Fig. ^| , 
and a fine structure provided tantalizing clues to the nature of such an effective 
string theory. A dramatic level reordering was observed as R became smaller, sug- 
gesting a bag model picture or multipole expansion may be more relevant at such 
scales. These energies were usecPI as a starting point in a Born-Oppenheimer treat- 
ment of heavy-quark mesons, both conventional and hybrid (bound by an excited 
gluon field). Agreement of level splittings from direct Monte Carlo calculations with 
those in the leading Born-Oppenheimer approximation validated such a treatment 
and provided a compelling physical picture of both conventional and hybrid heavy- 
quark mesons. 
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Fig. 1. The glueball mass spectrum-^! in the pure-gauge theory (left). The spectrum of gluonic 
excitationJSl in the presence of a static quark-antiquark pair separated by distance R (right). 



4. Baryons 

I am currently a member of the Lattice Hadron Physics Collaboration (LHPC). 
One of our goals is to compute the spectrum of baryon resonances, with an eye 
towards later determining the meson spectrum. We have designed l arge sets of ex- 
tended, gauge-invariant three-quark operators to facilitate this tastf^l. The usual 
operator construction which mimics the approach one would take in continuous 
space-time is very cumbersome, especially when tackling an entire spectrum. Our 
approach combines the physical characteristics of baryons with the symmetries of 
the lattice regularization of QCD used in simulations. For baryons at rest, our 
operators are formed using group-theoretical projections onto the irreducible repre- 
sentations (irreps) of the Oh symmetry group of a three-dimensional cubic lattice. 
There are four two-dimensional irreps G\ g , G% u , G2 g , G2 U and two four-dimensional 
representations H g and H u . The continuum- limit spins J of our states can be de- 
duced by examining degeneracy patterns across the different Oh irreps. 

® ®- •— ®— • I- *X* <T* 

single- singly- doubly- doubly- triply- triply- 

site displaced displaced-I displaced-L displaced- T displaced-0 

Fig. 2. The spatial arrangments of the extended three-quark baryon operators usecffSl. Smeared 
quark fields are shown by solid circles, line segments indicate gauge-covariant displacements, and 
each hollow circle indicates the location of a Levi-Civita color coupling. For simplicity, all displace- 
ments have the same length in an operator. 
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Fig. 3. Effective mass plots for three particular nucleon operators: a single-site operator in the Gi g 
channel (left), a doubly-displaccd-I operator in the Gi u channel (center), and a triply-displaced- 
T operator in the H g channel. For a correlation function C(t), the effective mass is defined by 
m effW = ln((7(<)/(7(t + 1)). These results are based on 25 quenched configurations on a 12 3 X 48 
anisotropic lattice using the Wilson action and optimized quark-field and link-variable smearing. 



Baryons are expected to be rather large objects, and hence, local operators will 
not suffice. Our approach to constructing extended operators is to use covariant 
displacements of the quark fields. Displacements in different directions are used to 
build up the appropriate orbital structure, and displacements of different lengths 
can build up the needed radial structure. There are six different spatial orientations 
that we use, shown in Fig [2 The singly-displaced operators are meant to mock up 
a diquark-quark coupling, and the doubly-displaced and triply-displaced operators 
are chosen since they favor the A-flux and Y"-flux configurations, respectively. 

We have now finished optimizing the quark-field and gauge-field smearing pa- 
rameters, and have begun low-statistics runs to prune out the ineffectual and overly- 
noisy operators. Some sample effective mass plots are shown in Fig. |3J Our goal is 
to find the smallest set of operators out of the several hundred we have constructed 
which are useful for extracting some number of lowest-lying states, and to push 
the technology to maximize the number of excited-state energies which can be re- 
liably determined. We shall report our findings in the near future. This work was 
supported by the U.S. National Science Foundation through grant PHY-0354982. 
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